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The production of quarkonium states below the open charm/bottom thresholds 
presents a particular challenge to QCD. Because of the relatively large quark masses, 
cc and hh production should be perturbatively calculable. However, the subsequent 
transition from the predominantly colour octet QQ pairs to physical quarkonium states 
can introduce nonperturbative aspects. These may lead to some model-dependence, 
requiring cross checks with as much data as possible. 

A generalisation of the colour evaporation model [1-5] provides a unihed approach 
to the production of the different quarkonium states below the open charm/bottom 
thresholds. As a specihc example, we consider charmonium production, although all 
arguments apply to bottonium production as well. Parton-parton interactions lead to 
the production of cc pairs, as shown in Fig. 1. We calculate the total “hidden” charm 
cross section, deg, by integrating over the cc pair mass from 2mc to 2m d- At high 
energy, the dominant production mechanism is gluon fusion (Fig. la), so that 

dcg(s) = / ds 

with g{x) denoting the gluon density and a the gg cc cross section. In pion-nucleon 
collisions, there are also signiheant quark-antiquark contributions (Fig. lb), which be¬ 
come dominant at low energies. Subsequently, the cc pair neutralizes its colour by 
interaction with the collision-induced colour held (“colour evaporation”). During this 
process, the c and the c either combine with light quarks to produce charmed mesons, 
or they bind with each other to form a charmonium state. More than half of the sub¬ 
threshold cross section deg in fact goes into open charm production (assuming rUc <1-5 
GeV); the additional energy needed to produce charmed hadrons is obtained (in general 
nonperturbatively) from the colour held in the interaction region. The yield of all char¬ 
monium states below the DD threshold is thus only a part of the total sub-threshold 
cross section: in this aspect the model we consider is a generalisation of the original 
colour evaporation model [1-5], which neglected the contribution of deg to open charm 
production. Using duality arguments, it equated deg to the sum over the charmonium 
states below the DD threshold. 

Neither the division of deg into open charm and charmonia nor the relative char¬ 
monium production rates are specihed by the generalised colour evaporation model. 
Hence its essential prediction is that the energy dependence of charmonium production 
is that of dcg(s). As a consequence, the ratios of different charmonium production cross 
sections are energy-independent. In Fig. 2, we show the ratio of J/'0 production from 
the decay Xc —^ 7 T/'0 to the total J/'0 production rate [5,6]. It provides a measure of 
the Xc/('^/'0) I'cife and is seen to be independent of incident energy for both pion and 
proton beams. In Fig. 3, we show the measured '0'/(J/'0) ratio [5,7,8]; it is also found 
to be independent of the incident energy, as well as of the projectile (pion or proton) 
and target (from protons to the heaviest nuclei [9]). Moreover, it is noteworthy that 
the ratio /{J/%1}) measured at high transverse momenta at the Tevatron [10] is quite 
compatible with the pr-integrated hxed target and ISR data (Fig. 4). 

The available bottonium data [11-14] also agree with constant production ratios, 
as seen in Fig. 5 for the ratios X’/T and T”/Tup to Tevatron energies. 


j dxidx 2 gixi) g{x 2 ) n(s) 5{s — x\X 2 s)^ (1) 
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The present data thus support one essential prediction of the colour evaporation 
model up to 1.8 TeV. We now check if it also correctly reproduces the variation of the 
production cross sections with incident energy in this region. In Figs. 6 and 7 we show 
the energy dependence of J/'0 production in pN collisions, 

^pJV^J/p('S) = fj/ip ^cc('S), (2) 

as obtained from the hidden charm cross section dec calculated in next-to-leading order 

[15] and with the normalisation hxed empirically. We have used the MRS D-’ 

[16] and GRV HO [17] parametrisations of the nucleon parton distributions functions 
[18]. For the GRV set, we have used rUc — 1.3 GeV, with both renormalisation and 
factorisation scales hxed to nic- In the MRS calculation, mc—1.2 GeV was used, with 
the scales set at 2mc- These parameters provide an adequate description of open charm 
production, although the results tend to lie somewhat below the measured total cc cross 
sections [19]. In Figs. 6 and 7 we show only the MRS D-’ result; the GRV HO result 
differs by less than 5 % in this energy range. The agreement with the data [5] over the 
entire range is quite satisfactory, with the normalisation = 0.025. In Fig. 8 we 
hud equally good agreement for the energy dependence of J/tf) production with pion 
beams. However, the fraction of J/'0 in the hidden charm cross section must be slightly 
higher to reproduce the pion data well, with = 0.034 for a good ht. This may 
well be due to greater uncertainties in the pionic parton distribution functions. We 
have also calculated the leading order cross section; the resulting theoretical K-factor, 
K = dcc^^/dc^^ remains between 2.0 and 2.5 over the currently measured energy range 
for both sets of parton distribution functions and for pion and proton beams. 

The fraction of dec producing charmonium rather than open charm is thus about 
10%. This is in accord with open charm calculations, which show [19] that much of the 
total cross section comes from subthreshold cc initial states which acquire the necessary 
energy for DD formation from the interaction colour held. To illustrate this. Fig. 9 
shows the fraction of the total open charm cross section with 2 nic < M <2 Mj:,^ where 
M = v^. It remains quite large even at very high incident energies. 

We further compare the longitudinal moment um dependence of charmonium pro¬ 
duction with recent experimental results. In hg. 10 we compare data with our calcula¬ 
tions for the xp dependence of J/^|^ production at several energies and for tt — p, p — p 
and p — p collisions [20]. Since there is a spread of integrated cross section values around 
the average deg, as seen in Figs. 6 - 8, we have normalised the calculated xp distribu¬ 
tion to the integrated experimental one. We conclude that the xp distributions are also 
consistent with the colour evaporation model. 

Next we comment briehy on the transverse momentum distributions. We are in¬ 
terested in low pt charmonium production, for which the model provides essentially 
no prediction. There is the intrinsic transverse momentum of the initial partons, the 
intrinsic momentum fluctuations of the colour held which neutralises the colour of the 
cc system in the evaporation process, and at larger pp higher order perturbative terms. 
Since there is no way to separate these different contributions in the low pp region, the 
model has no predictive power. 
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The colour evaporation model thus reproduces correctly both the energy depen¬ 
dence and the xp distributions of charmonium production, up to an open normalisation 
constant for each charmonium state, which can be hxed empirically by data. Once this 
is done, integrated and differential cross sections can be predicted for RHIC and LHC 
energies. From the hts to the data shown in Figs. 6 and 7 we obtain 


da. 


pN^J/tP 

dy 


y=0 


2.5 X 10“^ 


V dy 


( 3 ) 


for production. In Fig. 11 we show the resulting {dappi^jjp, jdy)y=Q as function 
of the center of mass energy, and in Fig. 12 we give the rapidity distributions at 
RHIC and LHC energies. The cross sections are listed in Table 1. 

Before commenting on our predictions, we hrst repeat the analysis for T production. 
Because the data generally give the sum of T, T’ and T” production, the branching 
ratios cannot simply be removed. Therefore we show in Fig. 13 the measured cross 
section for the sum of the three T states in the dilepton decay channel, denoted by 
B{da/dy)y=Q. We see that 


B 


da. 


pN^r 

dy 


y=0 


1.6 X 10“^ 


da^.^o 

00 

dy 


y=0 


( 4 ) 


gives a good description of the data up to and including ISR results. The results are 
also calculated using the MRS D-’ and GRV HO parton distribution functions, with 
mb = 4.75 GeV and the scales equal to mt- Assuming the bulk of the cross section to 
be from T(1S) production, and using the corresponding branching ratio, we estimate 
from eq. (4) that about 7% of the sub-threshold bb cross section leads to T production. 

Using the normalisaton determined in Eq. (4), we obtain the cross section for high 
energy T production; the results are shown in Figs. 14 and 15 and in Table 2. The recent 
high energy data from UAl [21] and CDF [14] agree very well with the predicted energy 
dependence, as seen in Fig. 13, giving strong support to the “new” parton distribution 
functions based on HERA data [22]. They also give us considerable co nfi dence in the 
extrapolation to LHC energy. 

We now comment on some features of our predictions. The two parton distribution 
functions, MRS D-’ and GRV HO, provide fully compatible results in the measured 
energy range. The T predictions agree with data even up to energies close to the LHC 
range. The MRS D-’ J/'ip cross section is about twice as large as the GRV HO prediction 
at LHC energies. This is because the MRS distributions require larger factorisation 
scales than the GRV distributions. Both parton distribution functions, with their chosen 
scales, also give acceptable hts to the measured open charm production cross sections 
(see [19] for more details). The difference thus gives some indication of the uncertainty 
of the J/'0 prediction. The Tresults agree over the entire energy range, since mb = 4.75 
GeV was used as the scale in both sets. 

At LHC energies, both the J/'ip and T rapidity distributions remain rather constant 
out to y ~ 4, using the MRS parton distributions. The GRV HO results show an even 
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wider plateau. In either case, there is a large window for forward detection at high 
energies. At RHIC energies, the J/'0 distributions are not as broad, with a forward 
plateau of 2 - 3 units for the MRS set, while the GRV distributions are somewhat 
narrower. The T rapidity distributions at RHIC energies are quite similar for both sets. 

Finally we note that the cross sections calculated with the recent parton distribution 
functions are considerably higher, typically by about a factor 20, than those given by 
an earlier empirical parametrisation, ~ exp(—15M/y^), labelled CR in Figs. 11 and 14 
[23-25] . This increase, conhrmed by new high energy T data (Fig. 14), is mainly due 
to the increase in the gluon distribution functions at small x, as suggested by data from 
HERA [22]. 

The colour evaporation model addresses the common energy behaviour of the dif¬ 
ferent quarkonium states. To determine their relative production rates, the colour evap¬ 
oration process has to be specihed in more detail. Let us consider one example of this. 
Assume that the initial colour octet state hrst neutralises its colour by interaction with 
the surrounding colour held, producing a colour singlet cc state. The relative weights of 
J/'0 and Ip' production can then be expressed [5] in terms of the corresponding masses 
and the squared charmonium wave functions at the origin, 

ff(</>') Rl'(o) / Y 
aW i?2(0) UV/ ' 

Here tp denotes the directly produced IS cc state, in contrast to the experimentally 
observed J/tp, 40% of which originates from radiative Xc decays (see Fig. 1). The 
wave functions at the origin can in turn be related to the dilepton decay widths Fee ~ 
(R2(0 )/m 2) [5], giving 

am 

am F^ ^ ^ 

Inserting the measured masses and decay widths, we hnd 


am 


~ 0.24 . 


(7) 


To compare this to the measured value of a{'ip')/a{'ip), we have to remove the Xe con¬ 
tributions from the experimental ratio. 


am) _ 

1 

■ am) 

am 

.7 ~ . 

_a{Jltp)_ 


( 8 ) 


With the experimental values o‘('0')/a( J/'0) ~ 0.14 (see Figs. 2 and 3) and {cr^^lajj^) ~ 
0.4 (see Fig. 1), this yields a('0')/cT('0) ~ 0.23, in good agreement with the theoretical 
result (7). We thus hnd that the projection of the colour singlet cc state onto J/^p and 
tp' correctly describes their production ratios at all energies and transverse momenta. 
The predictions for direct bottonium production ratios corresponding to eq. (7) are 


a(T) 


~ 0.36 


a{T") 


~ 0.27 . 


(9) 
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Since the contributions from indirect production through radiative Xb decay are not yet 
known and there is also feeding from higher S-states, a quantitative comparison is not 
possible here. Nevertheless, the predicted values differ from the data (see Fig. 5) by less 
than 50 % and hence appear reasonable. 

So far, the most complete description of the colour evaporation process is attempted 
in the colour singlet model [26], in which not only the cc formation but also the sub¬ 
sequent colour neutralisation is assumed to take place on a perturbative scale. The 
resonance formation is then determined by the appropriate wave functions with the 
right quantum numbers, as above. As a result, the production cross section for each 
charmonium state is completely determined to the order of perturbation theory used. 
Some characteristic production diagrams in lowest order are illustrated in Fig. 16. As 
generally formulated, the scale of the strong coupling constant in all perturbative dia¬ 
grams is determined by the mass of the heavy quark. 

Such a perturbative description of colour neutralisation can be justihed only if all 
moment um scales are sufficiently large. However, as seen in Fig. 16, colour neutralisa¬ 
tion for all but rjc, Xo c^nd and X 2 requires the emission or absorption of a “third” gluon. 
This restricts the possible applicability of the model to production at large transverse 
momentum. In the pT-integration, the “third” gluon is soft {k ~ Aqcd) in a signihcant 
part of phase space, and hence the model becomes unreliable here even though the in¬ 
tegration is infrared hnite. It is therefore not surprising that the colour singlet model 
leads to charmonium production ratios which disagree rather strongly with experiment. 
The quantum numbers of the X 2 allow partonic production at order ct^, while J/'0, Xi 
and 'Ip' production are all of order a^. Asa result, their production is much too strongly 
suppressed in comparison to the X 2 - Thus, while the model predicts X 2 l{Jj'>P) — 10, 
the measured ratios are below 2 [6]. The inclusion of certain relativistic corrections 
can somewhat reduce this discrepancy [5]. - Similar arguments hold for T produc¬ 
tion, although the soft part of the pT-iategration is relatively smaller, so that here the 
predictions may be closer to the data. 

For the validity of a perturbative treatment, the “third” gluon has to be hard 
enough to resolve the cc into individual quarks. Hence its momentum must be higher 
than 1/rj/^ ~ 1/(0.2 fm) ~ 1 GeV. Below this limit, it is not clear how colour neu¬ 
tralisation is achieved; presumably nonperturbative interactions of the colour octet cc 
with the gluon condensate play a considerable role. Some additional contributions can 
perhaps also be obtained by summing classes of perturbative contributions [27]. How¬ 
ever, as long as the additional interactions cannot be determined quantitatively, the 
prediction of the p^-integrated production ratios of the different charmonium states is 
not possible. It is not known if this is also true for bottonium production, or if here 
the role of soft processes has become sufficiently reduced to allow a fully perturbative 
treatment. 
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Table 1: J/^ Production 


^/s [GeV] 

{da/dy)f^^ [^b] 

(da/dy)^^^ [^b] 

20 

6.2 X 10-2 

5.8 X 10-2 


1.6 X 10-1 

1.5 X 10-1 

60 

2.5 X 10-1 

2.4 X 10-1 

100 

3.5 X 10-1 

3.4 X 10-1 

200 

6.3 X 10-1 

5.9 X 10-1 

500 

1.5 X lO+o 

1.2 X lO+o 

1000 

3.2 X lO+o 

2.5 X 10+° 

5500 

1.6 X 10+1 

5.9 X 10+° 

14000 

4.1 X 10+1 

1.1 X 10+1 


Table 2: (T + T' + T") Production 


y/s [GeV] 

{Bda/dy)f^^ [pb] 

(Bda/dy)^^^ [pb] 

15 

3.1 X 10-^ 

2.5 X 10-^ 

30 

9.7 X 10-1 

9.7 X 10-1 

60 

1.2 X 10+1 

1.2 X 10+1 

100 

3.4 X 10+1 

3.7 X 10+1 

200 

8.6 X 10+1 

1.0 X 10+2 

500 

2.5 X 10+2 

3.4 X 10+2 

1000 

5.5 X 10+2 

8.8 X 10+2 

5500 

3.0 X 10+3 

3.6 X 10+3 

14000 

7.8 X 10+3 

7.7 X 10+3 



















Figure Captions 


Fig. 1: Lowest order cc production through gluon fusion (a) and quark-antiquark anni¬ 
hilation (b). 

Fig. 2: The ratio of (xi + X 2 ) —^ J to total J/'0 production, as a function of the 
center of mass energy, ^/s^ by proton (open symbols) and pion beams (solid symbols) 
[ 6 ]. 

Fig. 3a: The ratio of '0' to J/%1) production as a function of the center of mass energy, 
^/s^ on proton (circles) and nuclear targets (squares) [5,7,8]. The average value is 
0.14 ± 0.03. 

Fig. 3b: The ratio of '0' to J/'0 production by proton beams as a function of the atomic 
mass number A for data in the energy range 20 < ^/s < 40 GeV [9]. The average value 
is 0.14 ± 0.01. 

Fig. 4a: The ratio of '0' to J/'0 production as a function of transverse momentum [10]; 
the shaded strip shows the average value of Fig. 3. 

Fig. 4b: The ratio of '0' to J/'0 production as a function of center of mass energy, ^/s. 
The hxed target and ISR data are integrated over the low px region, while the CDF 
point is the average over 5 < pt < 15 GeV. 

Fig. 5: The ratios of T’ and T” to T production as a function of the center of mass 
energy, ^/s [11-14]. The average values are 0.53T0.13 and 0.17T0.06, respectively. 

Fig. 6: The differential J/'ll) production cross section {da^JJ^/d'y) = 2.5 x 10“^ [da^^/dy) 
at p = 0, calculated with MRS D-’ parton distributions, compared to data [5]. 

Fig. 7: The J/'i/ production cross section = 2.5 x 10“^ iox xp > 0, calculated 
with MRS D-’ parton distributions, compared to data [5]. 

Fig. 8: The J/'ll) production cross section = 3.4 x 10“^ d'/A for xp > 0, calculated 
with MRS D-’/SMRS P2 parton distributions, compared to data [5]. 

Fig. 9: The fraction of the total open charm cross section due to the “hidden” charm 
mass interval [2mc,2m£)]. 

Fig. 10a: The J/'ll) longitudinal momentum distributions compared to pN and pN data 
[20], with xp = Pl{J/'/ 2)/pmaxiJ/'i/)] the MRS results are denoted by a solid, the GRV 
by a dashed line. 

Fig. 10b: The J/'i/ longitudinal momentum distributions compared to tiN data [20], 
with Xp = PL{ J/'>P)/Pmax{J/'>P)- the MRS results are denoted by a solid, the GRV by 
a dashed line. 

Fig. 11: Energy dependence of {da^J^/dy)y=o for J/'i/ production, as obtained with 
MRS D-’ and GRV HO parton distributions. 

Fig. 12: Rapidity distributions for J/'i/ production, calculated with MRS D-’ parton 
distributions at RHIG and LHG energies. 
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Fig. 13: The differential T production cross section {da^/dy) = 1.8 x 10“^ 

at y = 0, calculated with MRS D-’ parton distributions, compared to data [5]. The 

corresponding GRV HO predictions are very similar. 

Fig. 14: Energy dependence of {da^/dy)y=Q for T production, with high energy data 
from [21,14]; the predictions of MRS D-’ and GRV HO essentially coincide. Also shown 
(GR) is the phenomenological fit of [23]. 

Fig. 15: Rapidity distributions for T production calculated with MRS D-’ parton dis¬ 
tributions at RHIG and LHG energies. 

Fig. 16: Lowest order contributions to charmonium production in the colour singlet 
model. 
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